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RG. 4 Ran view of input buoyancy stnjcture and ftiutts (a), model outpu of 
f"1^Sir*°^'^ ^ topography and long- wavelength 

(> 1.000 km) geoid (c) for the model with vahaWe strikes of sJabs and 
surface velocrty and (d) for the model excluding the fautts. Contour 
interval for geoid is 10 m and dashed lines are negative values. The colour 
scales are the same as m Fig. 2. The length, width and thickness of the 
model boxes (in km) are 6.000. 7.500 and 1.500. respectively, which are 
covered by 54 x 54 x 16 elements. 



faults and slabs \vith variable strikes (Fig. 4^) shows that surface 
velocities remain plate-like (Fig. 4b) with P = 0.87 and 
Rt;p = 0.40 (Rj.p = 0.42 for a perfect plate). Surface velocities 
in R are predominantly parallel to the vertical faults and show 
oblique subduction where the dipping faults are not parallel to ihe 
spreading centre (Fig. 4^). Trench topography is evident ocean- 
ward from the dipping faults and there is 40 m long-wavclength 
geoid high over slabs (Fig. 4c). Qualitatively, both trenchyfore- 
buige topography and the geoid are essentially identical to 
observed features in the western Pacific (for example, approxi- 
mately 4 km deep and 150 km wide trenches, and 30tn long- 
wavelength geoid high). The reduced P results form the velocity 
variations m the spreading perpendicular direction near the short 
vertical fault (Fig. 4b) where slabs are younger; a more realistic 
temperature-dependent rheology may lead to a higher plalcncss. 
When the faults are excluded from the model, surface velocities 
vary gradually in both the spreading parallel and pcroendicular 
directions with F decreased to 0.75 (Fig. 4dl However. /It* only 
decreases slighUy to 0.39, suggesting that a large fraction of 
toroidal power is not a sufficient criterion for charaacriziM the 
existence of plates. 

These models show that the characteristics of plate motion can 
be explamed with realistic mantle buoyancy, weak faults and a 
power-law rheology (n = 3). Besides plate kinematics, the modeU 
also smiultancously predict realistic trench/fore-bulge topography 
and the long-wavelength geoid (Fig. 4b and c). In the past, 
different physical processes were used to explain these observed 
features (for example, elastic plate bending models for trench/ 
fore-bulge topography and laterally homogenous viscous flow 
models for the geoid). By simultaneously matching observed 
trench topography, plate kinematics and the long-wavelength 
geoid, our models will provide more realistic constraints on 
materials propenies and slab dynamics. Weak transform faults 
enhance the motion along the faults. Weak transform faults and 
their strong adjacent regions cause transform faults to guide plate 
motion (Fig. and h), consistent with inferences based on 
p ate kinematics'^ The guiding influence of transform faults on 
plate motion and the decoupling of thrust faults result in oblique 
subduction where the strike of thrust faults is non-perpendicular 
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to ^he transform faults (Figs 2h and 4*). consistent with observa- 
tions'-'*. Non-zero normal velocities on faults indicate migrating 
plate margins- Ultimately, with the formulation presented here 
and with multiple plates coupled with the solution of energy 
equation, we should be able to predict the evolution of plate 
geometry as accomplished in two dinlensions^ □ 
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ASTHMA now affects one child in seven in the United Kingdom*. 
Most cases (95%) of childhood asthma are associated with atopv, 
the immunogiobuJiD E (IgE) -mediated familial syndrome of 
allergic asthma, eczema and rhinitis. Segregation analysis has 
consistently suggested Ihe presence of major genes influencing 
atopy and IgE levels^, with the expectation that these genes may 
be identified by positional cloning or the examination of candi- 
date genes. Here we report the results of a genome-wide search for 
linkage to one qualitative and four quantitative traits associated 
with allergic (atopic) asthma. We have identified six potential 
Imkages {P < 0,001), five of which are to quantitative traits. 
Monte Carlo simulations show that 1.6 false-positive linkages 
at this level of significance would be expected from the data. One 
linkage, to chromosome llql3, has been established previously*. 
Three of the new loci show evidence of linkage to a second panel of 
families, in which maternal effects and pleiotropy of linked 
phcnotypes are seen. The results demonstrate the extent and 
the complexity of the genetic predisposition to asthma. 

Asthma is usually recognized epidemiologically by standard 
symptom questionnaires or by physician diagnosis*. The aon- 
specific bronchial hyper-responsiveness that accompanies 
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Marker 

D4S171 
DAS 1540 
D4S426 

D6S260 

06S276 

TNFA 

06S273 

D6S291 

D6S271 

D7S484 

D7S2250 

D7S528 

FCERB 
FGF3 
D11S96 
D11S901 

D13S262 
D13S270 
D13S153 
D13S170 

D16S421 
D16S515 
D16S516 
D16S289 
D16S507 
D16S505 



d 

0.01 
0.10 



0.21 
0.03 
0.03 
0.06 
0.10 



0.00 
0.03 



0.10 
0.05 
0.09 

0.04 
0.04 
0.20 

0.11 
0,12 
0.00 
0.06 
0.03 



Log^lgE 



P < 0.05 



asthma is measured by the dose 
response of airflow to a broncho- 
constrictor such as histamine or 
methachoiine. Atopy is detected 
by skin prick tests, by measure- 
ment of specific serum IgE titres 
against allergens using radio- 
allergosorbent test (RAST) or 
enzyme-linked immunosorbent 
assay (ELISA) techniques, or by 
quantifying the total serum IgE. 
Atopic disease is accompanied by 
blood eosinophilia, and eosinophils 

are prominent in asthmatic 

aiPAr'ays. 
The prevalence of atopy in 

Western populations is between 

40 and 50%'"^ The recruitment 

of families more than one 

asthmatic member results in 

samples in which 709t or more 

subjects are atopic, making it 

much more difficult to detect 

linkage'"'. For this reason, the 80 

families in the genome screen 

(from Busselton in Western 

Australia) were subselected 

from a population sample to 

include non-atopic members. 

We also reasoned that detection 

of linkage to quantitative traits 

would be enhanced by including 

subjects in normal and abnormal 

ranges of the trait distributions. 
We looked for linkage to four 

quantitative parameters: the total 

serum IgE. the skin test index . 
(STl), the peripheral blood eosinophil count, and bronchial 
responsiveness to methachoiine (slope)". A RAST index was 
also calculated, but gave no additional information to the STI 
and so was not included in the analysis described below. 

Differing indices of atopy may be increased in the same 
family'-", and RAST and skin-test responses reach a peak later 
in childhood than the total serum IgE, and decline at a slower rate 
thereafter'*. To account for this heterogeneity of phenotype 
(pleiotropy), the categorical trait of *atopy' was used, based on a 
combination of STI, RAST index, and total serum IgE. Linkage to 
quantitative traits was tested by the Haseman-Elston sib-pair 
technique'\ and to atopy by affected sib-pair methods. 

The 80 families sui\feyed contained a total of 203 offspring 
forming 172 sib-pairs. The mean age of the children was 12.6 ± 1.3 
(s.e.) years, their geometric mean IgE was 55.7 ± 1.1 kUr\ and 
their mean STI was 4.0 ± 0.41 mm; \2% of the children were 
asthmatic, and 25% admitted to wheeze; 52% were atopic, 36% 
were non-atopic, and 12% had an intermediate phenotype as 
defined; the latter were classified as unknown. The sib-pairs 
subdivided into 44 pairs in which both siblings were atopic. 87 
pairs in which one sibling was affected, and 41 in which both were 
unaffected or of intermediate phenotype. 

The families were initially screened with 269 markers (253 
autosomal and 16 X-linked). A genome map was constructed 
from these data, and the order and distances between the markers 
were compatible with previously published values'^ Based on 
these families, the markers had an average heterozygosity of 75%. 
and spanned approximately 3,100 centimorgans (cM). After initial 
statistical analysis, none of the X chromosome markers showed 
evidence of linkage to the five phenotype traits, and the remainder 
of this analysis therefore deals exclusively with autosomal maricers. 

Six regions of potential linkage (P < 0.001) to autosomal 
markers were detected with one or more phenotypes on chromo- 
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TABLE 1 Positive results of genome search (empirical P values) 



Skin test 
index 



P < 0.005 



P < 0 0005 P < 0.00005 

P < 0.005 
P < 0.01 



P < 0.05 
P < 0.05 
P < 0.05 
P < 0.0005 
P < 0.005 
P < 0.05 



Loge 
eosinophil 
count 



P < 0.05 
P < 0.0001 

P < 0.05 



P < 0.05 
P < 0.05 



Loge slope 

P < 0.05 
P < 0.05 
P < 0.0005 



Atopy 



P < 0.005 
P < 0.05 



P < 0.0005 
P < 0.05 
P < 0.05 



P< 0 05 
P< 0.005 
P< 0 001 
P < 0.05 



P < 0.05 

P < 0 05 
P < 0.05 



Only regions showing linkage to one or more phenotypes with P < 0.001 are shown. The recombination 
fraction id) is to the subsequent marker. 



somes 4, 6, 7, 1 1 . 1 3 and 16. respectively (Table 1 ). For the regions 
on chromosomes 4, 7 and 16. the test statistic met the more 
stringent criterion of P < 0.0005, and for two regions, on chromo- 
somes 6 and 1 1, P was <0.0001. Removal of sib-pair differences 
that were greater than 3 standard deviations from the mean did 
not significantly alter the P values for these loci. 

The preliminary identification of regions of potential linkage 
used correlated phenotypes with markers spanning most of the 
genome. Because of multiple testing, the true significance values 
will be greater than the normative values given in Table 1. To 
judge the global significance of the results, a Monte Carlo 
procedure was applied to calculate the number of regions that 
would be expected to reach these levels of significance in the 
absence of anv linkages (that is, the number of false-positive 
linkages). The procedure has advantages over traditional 
approaches to correcting nominal significance levels because it 
allows for the informativeness and distance between markers that 
were actually characterized in the genome-wide search. The 
expected number of regions of false-positive linkage was 1.6 for 
P < 0.001 (compared with 6 observed linkages), 0.8 for 
P< 0.0005 (5 observed linkages) and 0.2 for P < 0.0001 (2 
observed linkages), so the observed number of regions meeting 
these statistical criteria for linkage was much larger than would be 
expected by chance. In 35,000 simulations, P < 0.001 was 
observed in less than of replicates, indicating that the prob- 
ability of seeing six or more such regions without a true linkage was 
<0.005. 

Additional markers were studied on the six chromosomes that 
showed evidence for linkage (Table 1). Although the evidence in 
favour of linkage was not increased, in all instances at least two of 
the markers from each region had significance of P < 0.05. 

Chromosomes 11 and 16 exhibited potential linkage to IgE 
levels. On chromosome 1 1 the marker FCERB, a microsaiellite in 
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than exaci alleles were determined. Non-menderian mamer data were flagged by 
the programs GAS and UNKNOWN. Potentialty incorrect genotypes were re- 
examined and retested if necessary. Ts«o-point lod scores between nrtarKers and 
marker order were compared with published maps"-". 

Statistical analysis. Sib-pairs were included in the linkage analysis only when 
marker genotypes were known for both parents. Linkage was sought to atopy by 
affected s»b-pair analysis. Sibships containing more than one pair were 
weighted^. Linkage to quantitative trarts was calculated by the method of 
Haseman and Elston^*, m which the probabilities that sib-pairs have two, one 
or no alleles identical by descent are regressed o the square of the sib-pair trait 
differences^*. Linkage was examined further in the region of markers where P 
was <0.CX)1. For Table 1, linkage was presented wrthP < 0.01 or P < 0.05 if a 
marker is < 15 cM from a marker with a P < 0 001. The P values for linkage with 
the quantitative trait loci were assessed from the regression analysis of vanance 
(ANOVA), which assumes independence of the sib-pairs. In regions of potential 
linkage (P < 0.001). the P values were validated by empirical calculations in 
which thjs assumption is not required. The empirical statistics shown in Table 1 
were obtained by simulating the segregation of markers to offspring conditional 
on the observed parental genotypes without linkage to any of the traits. The 
empincal P values are the proportion of instances in which the linkage statistic 
equalled or exceeded its observed value in 500.000 simulated replicates. The 
fi^uency of false -positive linkages was assessed by a Monte Carlo simulation, 
in which replicate sets of observations were simulated under the absence of 
linkage with mer>delian segregation of the offspnng genotypes conditional on the 
observed parental genotypes. Linkage between markers was incorporated with 
recombination distances equal to those estimated from the data map obtained 
from this study, with offspring simulated under the assumption of lack of 
interference. The linkage statistics were calculated for each marker locus with 
the offspnng phenotype values as obsen^ed. The number of regions in which the 
linkage statistic for one or more phenotypes exceeded the assigned critical value 
was calculated over 35.000 replicates. Markers were assumed to detect 
different regions of linkage if all of the P values were non-significant over a 
region of >15cM between two markers with significant linkage. 

Replication. Panel B families were recruited from clinics in the United 
Klngdom^ Each family contained at least one asthmatic identified as a proband. 
Phenotyping was as for the Bussetton families, except that eosinophil counts had 
not been performed and bronchial responsiveness had not been measured. The 
families were tested with markers that showed potential linkage (P < 0.001) to 
any quantitative or qualitative trait in the Busselton families. Genotyping and 
linkage analysis were as for the Bussetton families. In addition, loci were tested 
for the presence of significant differences in linkage between maternally and 
paternally derived alleles: if a significam difference was found, the significance of 
linkage to maternal alleles was estimated. 
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AhfTERiOR- POSTERIOR (A-P) patterning is of fundamental impor- 
tance throughout vertebrate embr>'onic development Murine 
members of the triihorax (trx) and Potycomb group {Pc-Q) of 
genes regulate A-P patterning of segmented axial structures^'*, 
demonstrating conserved upstream regulation of homeotic path- 
ways between Drosophiia and mouse. Here we report the posi- 
tional cloning of a classical mouse mutation, eed (for embryonic 
ectoderm development)^ which is the highly conserved homologue of 
the Drosophiia Pc-G gene esc (for extra sex combs) ^ a long-term 
repressor of homeotic genes^ Mutants homozygous for a null 
allele of eed display disrupted A-P patterning of the primitive 
streak during gastrulation. Mutant embryos lack a node, noto- 
cfaord and somites, and there is no neural induction*. In contrast 
to absence of anterior structures, extra-embryonic mesoderm is 
abundant. Mice carrying a hypomorphic eed mutation exhibit 
posterior transformations along the axial skeleton. These results 
indicate that eed is required globally for A-P patterning through- 
out embryogenesis. 

The recessive, embryonic-lethal mutation eed maps to the 
albino-deletion complex on chromosome 7 (ref. 7), Previous 
genetic and molecular analyses positioned eed within a 150- 
kilobase (kb) intei^^al between the c^" and the c^^ deletion 
breakpoints' (Fig. \a). Here we have isolated eed from a 
CpG island (that is, a cluster of methylation-sensitive rare- 
cutter sites) located approximately 40 kb proximal to the c^^ 
deletion breakpoint''. Molecular access was obtained by 
chromosomal walking from the proximal end of yeast artificial 
chromosome (YAC) D6E7 (Fig. \a). A l.l-kbSflrt fragment from 
the CpG island was used to isolate overlapping clones from a 
mouse embryonic day 7.5 (E7.5) cDNA library generating a 1.95-kb 
cDNA (Fig. la). The candidate cDNA was validated as the eed 
gene by analysing two non-complementing ethyI-A/-nitrosurea 
(ENU)-induced alleles {l7Rn5'''''^ and /7/?n5'^*'^*) at the eed 
locus'^ The l7Rn5"''^^ and URnS'^"^'^ mutations represent 
null*^ and hypomorphic alleles'', respectively. DNA sequence 
analysis of reverse transcription- pol>'merase chain reaction 
(RT-PCR) products from homozygous E8.5 IJRnS'^^"*^^ embryos 
showed a transition of thymine at position 1040 to cytosine 
^1040 _^ ^pjg Yhis base exchange was absent in an 

E>nj-induced allele, l7Rn5'^^^^^, which /rarty-complements the 
l7Rn5^^^^^^ alIele^ Polymorphism can be excluded because both 
mutations were generated on the same BALB/cRl background^. 
As a second co-isogenic control, RT-PCR products from El 2.5 
homozygous l7Rn5^^^^ embryos also showed the expected wild- 
type T"***'. The sequence alteration in l7Rn5^'^^^ was identified as 
a T'°^^ — ^ A transversion which was absent in both l7Rn5^^^^^^ and 
l7Rn5^'''^^ (Fig. la). 

Conceptual translation revealed an open reading frame of 441 
amino acids, containing five tandemly repeated WD motifs (Fig. 
la). A potential PEST sequence" is defined by residues 22 to 52 of 
Eed protein (Fig. la). Two discrete clusters of basic amino acids 
encompassing a nine-residue spacer region (RK-X^-KKXKX, 
residues 65 to 80; Fig. la) are reminiscent of a bipartite nuclear 
targeting signal'\ Eed does not display any canonical DNA- 
binding motii. Three putative Asn-linked glycosylation sites 
(NXS/T)" at amino-acid positions 157, 283 and 349 (Fig. la) 
are compatible with a potential non-nuclear isoform of Eed 
generated by loss of the nuclear targeting signal through alter- 
native splicing (A. Schumacher and T. Magnuson, unpublished 
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the p-chain of the high-affinity receptor for IgE 
(FceRI-P), also showed linkage to the STI, but this 
was not observed on chromosome 16. The regions on 
chromosome 4 and 7 were linked to bronchial respon- 
siveness, whereas the region on chromosome 6 (near 
the class I genes of the major histocompatibility com- 
plex (MHC)) was linked to eosinophil counts. Weaker 
evidence of linkage {P < 0.01) was found with other 
phenotyp>es for the markers on chromosomes 6 and 7 
{Table 1). 

Markers from chromosome 13, around D13S153, 
showed evidence of linkage to the atopy phenotype in 
affected sib-pair analyses (Table 1). The sharing of 
maternal and paternal alleles identified by descent 
(IBD) was 46,1 and 20,0. 

The markers showing P < 0.001 for linkage were 
tested for replication in an additional group (panel B) 
of 77 nuclear and extended families, recruited from clinics in the 
United FGngdom, who had previously been used to map atopy on 
chromosome llql3^ These families had been characterized by 
the same protocols as the Busselton families, except that bronchial 
responsiveness and eosinophil counts had not been measured. 
They contained 215 offspring (268 sib-pairs), of which 61% were 
atopic and 56% asthmatic, the high values reflecting their origin 
through clinics. There were 121 atopic-affected sib-pairs, and 70 
asthmatic sib-pairs. 

Linkage of asthma to FCERB {P = 0.003) and to D16S289 
{P = 0.03) was seen in the panel B families (Table 2). Linkage 
with atopy was found to D13S153 (P = 0.003). The sib-pair 
sharing of Di3S153 alleles in these subjects was 65 IBD = 1, 
36 IBD = 0 (P = 0.003). The sharing for the combined data was 
1 1 1 IBD = 1, 56 IBD = 0. The combined Hodge-weighted statis- 
tic was 4. 158 (P = 0.000016). The empirical P value for linkage to 
D13SI53, calculated after 500,000 estimations, was <0.00001. 

Many studies have shown that maternal atopy and asthma are 
more likely to be transmitted to children than is paternal 
atopy*-"-^. We have not tested for maternal linkage in the 
genome screen families because of the difficulties induced by a 
further set of multiple comparisons. However, maternal linkage of 
atopy to chromosome llq markers has previously been demon- 
strated in the panel B families''^^ and so potentially linked 
markers identified in the genome screen were tested for parent 
of origin effects (Table 2). 

D4S426, FCERB and D16S289 showed significant differences 
in linkage between maternal and paternal alleles (Table 2). 
Linkage to maternal meioses was seen between D4S426 and 
atopy {P < 0.05) and the total serum IgE {P < 0.001), and 
between D16S289 and atopy {P < 0.01) and asthma 
(P < 0.001). As already observed, FCERB showed strong mater- 
nal linkage to atopy^^^ but a strong linkage to asthma was also 
seen (P < 0.00001) which has not been previously noted. The 
finding of maternal effects at several loci favours immunological 
interactions between mother and child, rather than genetic 
imprinting or anticipation, as a cause of these phenomena. 

It is not clear if linkage of different traits to different regions of 
the genome indicates that particular genes affect specific traits. 
However, a previously established region of genetic linkage, to 
chromosome llql3, may be used as an examplar*^ In the 
genome screen, the llql3 marker FCERB was linked to the IgE 
and the STL In other studies, FCERB or neighbouring markers 
have been linked to the composite phenotype of atopy'^ -'-^ and to 
bronchial responsiveness apparently independently of atopy". 
Similarly, D4S426 and D16S289 showed linkage to different 
phenotypes in genome screen and the panel B families. Although 
D13S153 showed linkage to the same phenotype of atopy in both 
Australian and British families, linkage between the total serum 
IgE and a polymorphism in the esterase D protein, near D13S153, 
has previously been suggested, with a lod score of 2.07 (ref. 25). 

The pleiotropy of phenotypes showing linkage to these loci are 
consistent with the known variability of atopy and asthma. Pleio- 
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TABLE 2 Testing for linkage in panel B families: combined and maternal meioses 



Locus IgE 


Skin test index Atopy 


Asthma 


D4S426 






(P < 0.001)* 


(P < 0.05) • 




D6S276 






D7S484 






FCERB 




P < 0.005 




(P < 0.0001)** 


(P< 0.00001) * 


D13S153 


P < 0.005 




D16S289 




P < 0.05 




(P < O.OD* 


(P < 0.001 )* 



Maternal meioses are given in parentheses. Statistical significance values are as 
follows: •P<0.01 patemal/matemal differences; ••P< 0.001 patemai/matemal 
differences. 



tropy may be attributable to environmental events, such as 
exposure to allergens, and to the age and selection of subjects. 
The relative distribution of quantitative traits in normal and 
abnormal ranges may also have differentially influenced the ability 
to detect linkage. 

The cytokine gene cluster on chromosome 5q31 has previously 
been reported to be linked to serum IgE in Amish families^*, and 
the result confirmed in Dutch asthmatic families-^ In our study we 
did not see linkage to either the serum IgE or to bronchial 
responsiveness with D5S393. D5S210 or D5S410 from the cyto- 
kine gene cluster, although we are currently typing multiple 
markers within this 20-cM region. 

In summary, both the Monte Carlo simulations and the replica- 
tion of some positive results in a second set of subjects suggest that 
we have identified regions of true linkage. Replication of linkage 
to complex traits is to be expected to be difficult*^, and testing for 
replication in other data sets and with a wider range of markers is 
still required for the linkages to chromosomes 6 and 7. It is 
unlikely that all important loci predisposing to atopy and 
asthma have been mapped, and larger studies and the pooling 
of data from several centres are required. □ 



Methods 

Subjects. The pnmary genome screen involved 364 subjects in 80 nuclear 
families subselected from a population sample of 230 families from Busselton 
in Western Australia''^. Clinical testing took place over three months in winter, 
to control for poiten exposure. Young families were recruited to reduce the effects 
of age and smoking on the lg£ levels and lung function. Families tn the genome 
screen irw;luded both atoptc and non-atopic members, and sibships of three or 
greater were not exclusively atopic or non-atopic. 

Phenotypes. Skin tests to house-dust mite (HDM) and mixed grass pollen 
(minus the response of negative controls), specific IgE litres to HDM and 
timothy grass, and the total serum IgE were measured^. A 'skin test index" (STl) 
was calculated as the sum of the skin-prick test results to HDM and grass mix; 
95% of individuals in this population who were atopic reacted to erther HDM or 
grass pollen, or both). Bronchial responsiveness to (up to 12 ^imol) methacho- 
line was measured^. The slope of the dose- response curve was calculated: 
(pre-dose forced expiratory volume in one second f FEVl) - last FEVl)/cumula- 
ttve dose of methacholine. A constant of 0.01 was added to each measurement, 
to allow log, transformation when the slope was <0. Eosinophils in peripheral 
blood swere Coulter-counted and the values log,-transformed before analysis. 
Atopy was defined as STl > 5 mm, or a RAST score to HDM and timothy grass 
>2, or a total senjm IgE > the 7th decile of the age-corrected population. 
Normal' was defined as an STl of 0 and a RAST index of 0, and a total serum 
IgE < the 7th decrle of the age- and sex-matched population. Intermediate 
phenotypes were classified as unknown. The subjects were administered a 
modified Brrttsh MRC questionnaire^®. Asthma was defined as a positive answer 
to the questions '*Have you ever had an attack of asthma?" and 'If yes. has this 
happened on more than one occasion?" 

QefKrtypes. Genomic ONA was isolated from penpherai blood leucocytes by 
phenoVchtoroform extraction. We typed 22 dtnculeotide repeat mart<ers by 
radioactive techniques^* (marker names italicized in Table 1). For 274 markers. 
50 ng of each DNA sample was genotyped by fluorescence-based semi-auto- 
matic methods Allele sizes were determined using the GENOTYPER soft- 
wares^. Data were converted into numbered alleles using the GAS (version 2.0) 
program (A. Young, unpublished). Allele size differences within families rather 
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